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We have measured the polarizations of J/ip and tp(2S) mesons as functions of their transverse 
momentum pr when they are produced promptly in the rapidity range \y\ < 0.6 with pr > 5 GeV/c. 
The analysis is performed using a data sample with an integrated luminosity of about 800 pb _1 
collected by the CDF II detector. For both vector mesons, we find that the polarizations become 
increasingly longitudinal as pr increases from 5 to 30 GeV/c. These results are compared to the pre- 
dictions of nonrelativistic quantum chromodynamics and other contemporary models. The effective 
polarizations of J/yj and yj(2S) mesons from i?-hadron decays are also reported. 
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An effective field theory, nonrelativistic quantum chrome-dynamics (NRQCD) provides a rigorous formalism for 
calculating the production rates of charmonium (cc) states. NRQCD explains the direct production cross sections 
for J/ip and ip(2S) mesons observed at the Tevatron 0, 0] and predicts their increasingly transverse polarizations 
as pt increases, where px is the meson's momentum component perpendicular to the colliding beam direction 0]. 
The first polarization measurements at the Tevatron || did not show such a trend. This Letter reports on J/ip and 
ip(2S) polarization measurements with a larger data sample than previously available. This allows the extension of 
the measurement to a higher px region and makes a more stringent test of the NRQCD prediction. 

The NRQCD cross section calculation for cc production separates the long-distance nonperturbative contributions 
from the short-distance perturbative behavior. The former is treated as an expansion of the matrix elements in powers 
of the nonrelativistic charm-quark velocity. This expansion can be computed by lattice simulations, but currently the 
expansion coefficients are treated as universal parameters, which are adjusted to match the cross section measurements 
at the Tevatron 0, The calculation also applies to cc production in ep collisions, but HERA measurements of J/ip 
polarization tend to disagree with the NRQCD prediction Q. These difficulties have led some authors to explore 
alternative power expansions of the long-distance interactions for the cc system Q. There are also new QCD-inspired 
models, the gluon tower model [§| and the fcr-factorization model ||, that accomodate vector- meson cross sections at 
both HERA and the Tevatron and predict the vector- meson polarizations as functions of px- These authors emphasize 
that measuring the vector-meson polarizations as functions of px is a crucial test of NRQCD. 

The CDF II detector is described in detail elsewhere In this analysis, the essential features are a muon system 

covering the central region of pseudorapidity, \rj\ < 0.6, and the tracking system, immersed in the 1.4 T solenoidal 
magnetic field and composed of a silicon microstrip detector and a cylindrical drift chamber called the central outer 
tracker (COT). The data used here correspond to an integrated luminosity of about 800 pb _1 and were recorded 
between June 2004 and February 2006 by a dimuon trigger, which requires two opposite-charge muon candidates, 
each having pt > 1.5 GeV/c. 

Decays of vector mesons V (either J/ip or ip(2S)) — » are selected from dimuon events for which each track 

has segments reconstructed in both the COT and the silicon microstrip detector. The px of each muon is required 
to exceed 1.75 GeV/c in order to guarantee a well-measured trigger efficiency. The muon track pair is required to 
be consistent with originating from a common vertex and to have an invariant mass M within the range 2.8 (3.4) < 
M < 3.4 (3.9) GeV/c 2 to be considered as a J/ip (ip(2S)) candidate. To have a reasonable polarization sensitivity, 
the vector-meson candidates are required to have pt > 5 GeV/c in the rapidity range \y (= ^ln ^^|| )| < 0.6, where 
E is the energy and pn is the momentum parallel to the beam direction of the dimuon system. Events are separated 
into a signal region and sideband regions, as indicated in Fig. [T] The fit to the data uses a double (single) Gaussian 
for the J/ip (ip(2S)) signal and a linear background shape. The fits are used only to define signal and background 
regions. The signal regions are within 3oy of the fitted mass peaks My, where ay is the width obtained in the fit 
to the invariant mass distribution. Both the background distribution and the quantity of background events under 
the signal peak are estimated by events from the lower and upper mass sidebands. The sideband regions are 7<Jj/^, 
(4cr^(2s)) away from the signal region for J/ip (ip(2S)). 

For each candidate, we compute ct = M L xy /px, where t is the proper decay time and L xy is the transverse distance 
between the beam line and the decay vertex in the plane normal to the beam direction. The ct distributions of the 
selected dimuon events are shown in Fig. O The ct distribution of prompt events is a Gaussian distribution centered 
at zero due to finite tracking resolution. For J/ip, the prompt events are due to direct production or the decays of 
heavier charmonium states such as Xc an d ip(2S); for ip(2S), the prompt events are almost entirely due to direct 
production since heavier charmonium states rarely decay to tp(2S) [llj. Both the J/ip and the ip(2S) samples contain 
significant numbers of events originating from long-lived B-hadron decays, as can be seen from the event excess at 
positive ct. We have measured the fraction of B — > J/ip + X events in the J / ip sample and found agreement with 
other results [3| . We select prompt events by requiring the sum of the squared impact parameter significances of the 

positively and negatively charged muon tracks S = {^-) 2 + (^-) 2 < 8. The impact parameter do is the distance 
of closest approach of the track to the beam line in the transverse plane. Vector-meson candidates from B-hadron 
decays are selected by requiring S > 16 and ct > 0.03 cm. This requirement retains a negligible fraction of prompt 
events in the B sample. 

To measure the polarizations of prompt J/ip and ip(2S) mesons as functions of px, the J/ip events are analyzed in 
six pt bins and the ip(2S) events in three bins, shown in Table [T] We determine the fraction of -B-decay background 
remaining in prompt samples fbkd by subtracting the number of negative ct events from the number of positive ct 
events. Only a negligible fraction (< 0.2%) of B decays produce vector-meson events with negative ct. For both 
vector mesons, fbkd increases with px, as listed in Table |U The prompt polarization from the fitting algorithm is 
corrected for this contamination. 
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FIG. 1: Invariant mass distributions for (a) J/tp and (b) ip(2S) candidates. The curves are fits to the data. The solid (dashed) 
lines indicate the signal (sideband) regions. 
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FIG. 2: Sideband-subtracted ct distributions for (a) J/tp and (b) ip(2S) events. The prompt Gaussian peak, positive excess 
from _B-hadron decays, and negative tail from mismeasured events are shown. The dotted line is the reflection of the negative 
ct histogram about zero. 



The polarization information is contained in the distribution of the muon decay angle 9* , the angle of the /i + in the 
rest frame of vector meson with respect to the vector-meson boost direction in the laboratory system. The decay angle 

d N 

distribution depends on the polarization parameter a: oc 1 + a cos 2 9* (—1 < a < 1). For fully transverse 

a cos 6 

(longitudinal) polarization, a = +1 (—1)- Intermediate values of a indicate a mixture of transverse and longitudinal 
polarization. 

A template method is used to account for acceptance and efficiency. Two sets of cos 6* distributions for fully 
polarized decays of J/tp and ij>(2S) events, one longitudinal (L) and the other transverse (T), are produced with the 
CDF simulation program using the efficiency-corrected pt spectra measured from data [3|, 1 1 21 ] . We use the muon 
trigger efficiency measured using data as a function of track parameters (pr, rj, 4>) to account for detector non- 
uniformities. The parametrized efficiency is used as a filter on all simulated muons. Events that pass reconstruction 
represent the behavior of fully polarized vector-meson decays in the detector. 

The fitting algorithm Q uses two binned cos 6* distributions for each px bin, one made by N$ events from the 
signal region (signal plus background) and the other made by Ng events from the sideband regions (background). 
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The \ 2 minimization is done simultaneously for both cos 9* distributions. The fitting algorithm includes an individual 
background term for each cos 9* bin, normalized to Nb ■ Simulation shows that the cos 9* resolution at all decay 
angles over the entire px range is much smaller than the bin width of 0.05 (0.10 for xp(2S)) used here. The data, fit, 
and template distributions for the worst fit (9% probability) in the J/ip data are shown in Fig.[3J 

7 < p T (GeV/c) < 9 

0.05 



0.04 
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"-1 -0.8 -0.6 -0.4 -0.2 ,0.2 0.4 0.6 0.8 1 
COS 6 

FIG. 3: cos 9* distribution of data (points) and polarization fit for the worst x 2 probability bin in the J/ip data. The dotted 
(dashed) line is the template for fully L (T) polarization. The fit describes the overall trend of the data well. 

All systematic uncertainties are much smaller than the statistical uncertainties. Varying the pr spectrum used in 
the simulation by la changed the polarization parameter for J/ip at most by 0.002. A systematic uncertainty of 0.007 
was estimated by the change in the polarization parameter when a modification was made on all trigger efficiencies 
by ±1<7. For ip(2S), the dominant systematic uncertainty came from the yield estimate because of the radiative tail 
and the large background. The total systematic uncertainties shown in Table Q] were taken to be the quadrature 
sum of these individual uncertainties. Other possible sources of systematic uncertainties - signal definition and cos 9* 
binning - were determined to be negligible. Corrections to prompt polarization from i?-decay contamination were 
small, so that uncertainties on £?-decay polarization measurements also had negligible effect. No (^-dependence of the 
polarizations was observed. 

The polarization of J/ip mesons from inclusive B u and Bd decays was measured by the BABAR collaboration 
In this analysis, the -B-hadron direction is unknown, so we define 9* with respect to the J/ip direction in the laboratory 
system. The resulting polarization is somewhat diluted. As discussed in Ref. CDF uses a Monte Carlo procedure to 
adapt the BABAR measurement to predict the effective J/ip polarization parameter. For the J/tp events with 5 < pr < 
30 GeV/c, the CDF model for B u and Bd decays gives a e jf = —0.145 ± 0.009, independent of pr- We have measured 
the polarization of vector mesons from B-hadron decays. For J/ip, we find <x e ff = — 0.106±0.033 (stat)± 0.007 (syst). 
At this level of accuracy, a polarization contribution by J/tp mesons from B s and 6-baryon decays cannot be separated 
from the effective polarization due to those from B u and Bd decays. We also report the first measurement of the 
i/j(2S) polarization from B-hadron decays: a e ff = 0.36 ± 0.25 (stat) ± 0.03 (syst). 

The polarization parameters for both prompt vector mesons corrected for f^kd using our experimental results on a e f / 
are listed as functions of pr in Table Q] and are plotted in Fig. [4] The polarization parameters for J / ip are negative 
over the entire pr range of measurement and become increasingly negative (favoring longitudinal polarization) as 
Pt increases. For ip(2S), the central value of the polarization parameter is positive at small pr, but, given the 
uncertainties, its behavior is consistent with the trend shown in the measurement of the J/ip polarization. 

The polarization behavior measured previously with 110 pb _1 Q is not consistent with the results presented here. 
This is a differential measurement, and the muon efficiencies in this analysis are true dimuon efficiencies. In Ref. [f|, 
they are the product of independent single muon efficiencies. The efficiency for muons with px < 4 GeV/c is crucial for 
good polarization sensitivity. In this analysis, the muon efficiency varies smoothly from 99% to 97% over this range. 
In the analysis of Ref. it varied from 93% to 40% with significant jumps between individual data points. Data 
from periods of drift chamber aging were omitted from this analysis because the polarization results were inconsistent 
with the remainder of the data. Studies such as this were not done in the analysis of Ref. [f| . The systematics of the 
polarization measurement are much better understood in this analysis. 

These polarization measurements for the charmed vector mesons extend to a pr regime where perturbative QCD 
should be applicable. The results are compared to the predictions of NRQCD and the fey-factorization model in Fig. |U 
The prediction of the fey-factorization model is presented for pt < 20 GeV/c and does not include the contribution 
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p T (GeV/c) <p T > (GeV/c) 


fbkd{%) 






5- 


-6 


5.5 


2.8 ±0.2 


-0.004 ± 0.029 ±0.009 15.5/21 




6- 


-7 


G.5 


3.4 ±0.2 


-0.015 ± 0.028 ± 0.010 24.1/23 




7- 


-9 


7.8 


4.1 ±0.2 


-0.077 ± 0.023 ± 0.013 35.1/25 




9- 


12 


10.1 


5.7 ±0.3 


-0.094 ± 0.028 ±0.007 34.0/29 




12- 


-17 


13.7 


6.7 ±0.6 


-0.140 ±0.043 ± 0.007 35.0/31 




17- 


-30 


20.0 


13.6 ± 1.4 


-0.187 ± 0.090 ± 0.007 33.9/35 




5- 


-7 


5.9 


1.6 ±0.9 


+0.314 ± 0.242 ± 0.028 13.1/11 




7- 


10 


8.2 


4.9 ± 1.2 


-0.013 ±0.201 ± 0.035 18.5/13 




10- 


-30 


12.6 


8.6 ± 1.8 


-0.374 ± 0.222 ± 0.062 26.9/17 



TABLE I: Polarization parameter a for prompt production in each pr bin. The first (second) uncertainty is statistical (sys- 
tematic). <pr> is the average transverse momentum. 
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FIG. 4: Prompt polarizations as functions of pr- (a) J/ip and (b) ip(2S). The band (line) is the prediction from NRQCD Q 
(the fcT-factorization model ||). 



from the decays of heavier charmonium states for J/ip production. The polarizations for prompt production of 
both vector mesons become increasingly longitudinal as pr increases beyond 10 GeV/c. This behavior is in strong 
disagreement with the NRQCD prediction of large transverse polarization at high px- It is striking that the NRQCD 
calculation and the other models reproduce the measured J/ip and ip{2S) cross sections at the Tevatron, but fail to 
describe the polarization at high px- This indicates that there is some important aspect of the production mechanism 
that is not yet understood. 
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